Bose-Einstein condensation in antiferromagnets close to the saturation field 
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At zero temperature and strong applied magnetic fields the ground sate of an anisotropic antiferro- 
magnet is a saturated paramagnet with fully aligned spins. We study the quantum phase transition 
as the field is reduced below an upper critical H C 2 and the system enters a XY-antiferromagnetic 
phase. Using a bond operator representation we consider a model spin-1 Heisenberg antiferromag- 
netic with single-ion anisotropy in hyper-cubic lattices under strong magnetic fields. We show that 
the transition at H C 2 can be interpreted as a Bose-Einstein condensation (BEC) of magnons. The 
theoretical results are used to analyze our magnetization versus field data in the organic compound 
NiCh-4:SC(N £[2)2 (DTN) at very low temperatures. This is the ideal BEC system to study this 
transition since H C 2 is sufficiently low to be reached with static magnetic fields (as opposed to 
pulsed fields) . The scaling of the magnetization as a function of field and temperature close to H C 2 
shows excellent agreement with the theoretical predictions. It allows to obtain the quantum critical 
exponents and confirm the BEC nature of the transition at H C 2- 
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The organic compound NiCl 2 -±SC{N H 2 ) 2 (DTN) 
undergoes a field induced non-magnetic to XY- 
antiferromagnetic transition 1 ^. This transition can be 
viewed as a Bose-Einstein condensation of magnons as- 
sociated with the Ni spin 1 degrees of freedom. Other 
magnetic systems with a singlet ground state either with 
spin-1 Ni atoms or spin- 1/2 dimers have also been shown 
to exhibit this transitions^. At zero temperature it is 
driven by the magnetic field H that reduces the Zeeman 
energy of the S z — 1 state until it becomes degenerate 
with that of the product state of S* = 0. At this point, 
H = H c i, the antiferromagnetic (AF) interactions give 
rise to a long range ordered phase. Experimentally, the 
magnetization M at very low temperatures starts to in- 
crease above the critical magnetic field H c \ and eventu- 
ally saturates above a critical magnetic field fl C 2«i£iL2j£. 
The transition at H cl has been intensively investigated, 
both theoretically 1 --^ and experimentallyi^^ 7 -^, while 
the one at H c2 is much less studied. The DTN is the ideal 
BEC system to investigate the latter transition since de- 
tailed magnetization curves can be obtained close to the 
critical field H c2 = 12. 3T. In other well known BEC sys- 
tems, as BaCuSi 2 Ogl2- and TlCuCl^, the critical fields 
H c2 are 49T and 83T, respectively, and presently can only 
be reached using pulsed fields. The excellent quality of 
the magnetization versus field curves obtained in DTN 
using standard superconducting coils is essential for the 
scaling analysis presented here. 

In this paper we study the transition at H c2 . A theo- 
retical approach is more directly developed, starting from 
the saturated paramagnetic (PARA) phase. We consider 



decreasing the external magnetic field at T = to the 
critical value H c2 where the transverse components of 
the magnetization condense. A scaling approach for this 
transition has recently been proposed^. In this Com- 
munication we provide the microscopic theory for this 
transition. We identify its universality class as a Bose- 
Einstein condensation associated with a dynamic expo- 
nent z = 2. We compare the predictions for the scal- 
ing behavior of the magnetization close to the quantum 
critical point (QCP) (T = 0, H = H c2 ) with experimen- 
tal magnetization data on DTN and obtain an excellent 
agreement. 

For a long time the magnetically ordered state and 
low energy excitations of quantum Heisenberg magnets 
have been studied using the spin-wave expansion (see^ 
and references therein). This is usually implemented by 
expressing the components of the spin operators at lat- 
tice sites i in terms of canonical boson operators b k and 
bk using the Holstein-Primakoff (H-P) transformation 16 , 
the Dyson-Maleev (D-M) transformatio n 17 ' 18 , or the 
Schwinger transformation 19 (ST). Here we use the bond- 
operator mean-field theory (BOMFT) 20 to study spin- 
1 Heisenberg AF in hyper-cubic lattices with single-ion 
anisotropy close to the quantum phase transition at the 
saturation field H c2 . It yields the phase diagram and the 
thermodynamic behavior of the model close to H c2 . The 
BOMFT gives an exact description of this transition for 
three dimensional (3d) systems. The reason is that the 
effective dimension d e f / = d + z = 5 associated with the 
QCP is larger than the upper critical dimension d c = 4 
above which mean-field theory is exact. 
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The Hamiltonian describing the magnetic system is, 
H=±J2( S ? S * + S t S j+ S * S j) + DJ2(S*) 2 -HJ2 S? 

(1) 

where the sum is over all nearest neighbor pairs of a d- 
dimensional hyper-cubic lattice with N sites occupied by 
spins with S = 1. J > is the AF exchange coupling, D 
is the single-ion anisotropy, H the magnetic field applied 
in the z-direction (gfi B — 1)- Starting from the bond- 
operator representation for two spins S = 1/2 20 , Wang 
and collaborator o 12 ' 13 i 21 obtained a representation for a 
spin-1 Heisenberg system with a single- ion anisotropy in 
terms of these operators. At zero temperature and for 
external magnetic fields larger than the saturation mag- 
netic field H C 2 the spins are fully aligned with the field. 
In this case the bond operator representation can be ex- 
pressed as, 



S+ = V2ut Zl S~ = V2ut\, S z = 1 



(2) 



with the constraint u\ui 



t f t 



1. We have used 
that in this large field case the components of the spins 
perpendicular to the field can be projected out and those 



parallel condense, such that, Ui 



,t _ 



u. This mapping 



is exact for H > H C 2, where the probability of the down 
spin state d'd\0) is strictly zero at T = 0. The mag- 
netic ordering for H = H C 2 can be identified as a Bose- 
Einstein condensation of the transverse components of 
the spins which give rise to the collective magnon ex- 
citations. The operators t\ z U tZ describe the departure 
of the spins from the field direction and are associated 
with these excitations. Replacing Eq. ([2]) in Eq. (fl]) 
with S 1 * = S x ± iS v and changing from atomic to nor- 



mal coordinates, i.e., taking t\ z — -^J2k' 
have, 



ik.ri lA 



= ^^Ukb\b k + E g , 



bl, we 



(3) 



where, 



Eg = N + Du 2 -H- {in 2 + f?j (4) 

is the ground state energy of the system. The dispersion 
relation of the excitations for H > H C 2 is given by, 



H -(£> + fj,)- JZ{l-u 2 lk ) 



(5) 



with, 7^ = d^ 1 Y^t=\ cos(k.a„) and Z the number of 
nearest neighbors. The chemical potential /i was intro- 
duced to impose the constraint condition of single oc- 
cupancy. This and the parameter u are determined by 
solving coupled, self-consistent, saddle point equations^ 3 -. 

For fixed field H or temperature T, the thermody- 
namic quantities can be obtained from the internal en- 
ergy U. This is given by^ 3 -: U = J2k < b\b k > +E g , 



with the Bose factor, n k =< b\b k >= | ^coth — lj 

and E g given by Eq. (U). For simplicity, we impose 
boundary conditions on a d-dimensional hyper-cubic lat- 
tice with primitive lattice vectors and lattice spac- 
ing a — |a„| = 1. The minus sign in front of J 
takes into account explicitly the antiferromagnetic na- 
ture of the exchange interactions. For hyper-cubic lat- 
tices, the minimum of the spin-wave spectrum occurs 
for q = Q = (ir/a, it /a, n/a) in three dimensions. Since 
ujq = H — (D + fx) — JZ(u 2 + 1), the condition loq = 
defines the critical field H c2 = D + /i + JZ(u 2 + 1), below 
which the spin-wave energy becomes negative signaling 
the entrance of the system in the AF phase. Finally, writ- 
ing k = Q + q and expanding for small q, the spin- wave 
dispersion relation can be obtained as follows: 



uj q = (H — H c2 ) + Dq 2 , 



(6) 



where the spin-wave stiffness at T = is given by, 
D = Ju 2 . The quantum phase transition at H C 2 has 
a dynamic exponent z = 2 due to the ferromagnetic-like 
dispersion of the magnons, in spite of the antiferromag- 
netic character of the interactions^. 

Magnetization — Close to the critical field H C 2 the 
temperature-dependent magnetization should follow a 
power law^. We define the variation of uniform mag- 
netization per volume V as, AM = (M sat — M)/V = 
Sfe < ^\j )k > where M sat is the saturation magnetiza- 
tion. Considering the spectrum of excitations, Eq. ((5]), 
we have in the thermodynamic limit, 



AM: 



Sd 



(k B Tf/ 2 dx(x-yf- l (coth^-l) 



(7) 

where x = I3u q = y + (3 2 T)q 2 , y — (38, Sd the solid angle 
and we have defined 5 =\H — H C 2\ as the distance to the 
QCP. At this point, we have to consider in which region 
of the phase diagram (see Fig.l of RefJ^) we are inter- 
ested. Because we want to calculate the magnetization 
above H C 2 and particularly at the quantum critical tra- 
jectory, H = H C 2, T — > 0, we consider region II in Figure 
1 of Ref.— where k B T ^> 8 and consequently y <C I . 
Calculating the integral above in three dimensions, we 
obtain, AM 3d =polylog(f , e"") (k B T)^ 2 / (n 3 / 2 ® 3 / 2 ) , 
where polylog(a, z) = Y^Li z ™ / n ° ^ s the general poly- 
logarithm function of index a at the point z. Along the 
quantum critical trajectory, H = H C 2, T — > 0, we find 
AM 3d (S - 0) - ^(3/2) (fc s T) 3 / 2 / {it 3 ' 2 ® 3 ' 2 ), where 
C, is the Riemann zeta-function. 

Specific heat — From the internal energy U obtained 
before and using the thermodynamic relation, Cy = 
OU/dT we get, 



k B S d (k B T) d / 2 
2n d D d / 2 



dxx 2 (x — y) 



F_1 sinh" 2 



(I)- 



. . (8) 

Again, we consider y -C 1 and along the quantum critical 
trajectory in 3d, we find, CV(<5 = 0) = 15fcs(-\/i -I- An — 
\){k B Tf 2 /(8tt 3 / 2 S) 3 / 2 ). 



Susceptibility — Here we use the relation, \ — dM/dH, 
where M is the magnetization. Taking the derivative of 
Eq. ((7j) with respect to H and changing variables, we 
have at S = 0, 



X 



s d 



-{k B Tf. 



dxx 2 1 sinh 



(!)■ < 9 » 



For 3d the longitudinal susceptibility at the critical field 
is given by, X (<5 = 0) = (fc B T) 1 / 2 /((vr - l)^® 3/2 )- 

The temperature dependence of the physical proper- 
ties calculated above in 3d, can be easily obtained from 
the fact that the QCP at H c2 is governed by Gausssian 
exponents and the free energy has a scaling form 24 , 



f<x\6\ 2 - a F(T/\6n, 



(10) 



where 8 = (H — H c2 ) as defined before. The Gaussian 
nature of the exponents in 3d is a consequence that the 
effective dimension d e ff = d + z = 5, which is larger than 
the upper critical dimension^, d c = 4. 

Critical line — The critical line that separates the po- 
larized PARA state from the AF phase with nonzero 
staggered magnetization can be written on the neighbor- 
hood of the QCP as, T N (H) cx \H c2 - . Theories 
for a 3d Bose ga o 5 ' 25 and mean-field treatment 6 - give a 
universal value, ip = 2/3. Scaling theory shows that 
although d e ff > d in three dimensions, the magnon- 
magnon interaction is dangerously irrelevant and must 
be considere d 24 ' 26 . Then we expect that the quartic cor- 
rections to the mean-field result Eq. ([3]) will be now 
important. Within the H-P representation for the spin 
operators, the mean-field Hamiltonian including the dy- 
namical spin wave interactions is given by, 



k,k' 



JZ 



D 



bUkb 



E„ 



where (4p + D) works as an effective repulsion between 
the magnons. The dispersion Uk is given by Eq. (O and 
Eg by Eq. (fj|. We decouple the spin- wave interaction 

as, b\bkb\.,bk> ~ b\bk < b\bk >■ Thus, trivially we obtain 
the internal energy as, 



^u/ q < b\b q > 



-E, 



u ■ 



(12) 



where we have already considered the proximity to the 
QCP (H = H c2 , T = 0). The spectrum of excitations 
taking into account magnon-magnon interactions is, 



(H-H c 



JZ 



D < blbk > +3V- (13) 



We set up an equation for the critical line within the 
mean-field approximation where the effect of magnon- 
magnon fluctuations are included in a self-consistent 
manner. The critical temperature Tn(H) is determined 
by the condition 6(H,Tn) — where, 

S(H,T N ) = {H-H c2 ) + 



JZ_ 



-D 



Sd 
4T 2 



d/2 



dxx^ 1 (coth 1 



(14) 

For 2d the integral above diverges as expected from 
general arguments 2 ^. For 3d we get, IcbTn — 
(C(3/2)(JZ/2 + D))- 2 / 3 ttD(H c2 -H) 2 / 3 . Notice that the 
effective magnon-magnon coupling strength (JZ/2 + D) 
determines the transition temperature, despite the Gaus- 
sian exponents, as expected from the dangerously irrel- 
evant nature of the magnon-magnon interactions. If we 
write the equation for the critical line, 5(H,T) — 0, in 
the form, H c2 (T) = H c2 (0) - vqT 1 ^, with v related 
to the spin-wave interaction, we identify the shift expo- 
nent, ip = zj(d + z — 2) = 2/3, in agreement with the 
renormalization group (RG) result^ 6 .. The temperature 
dependence of 6 arising from the spin-wave interactions 
can modify the temperature dependence of AM, Cy and 
X at H = H c2 . In the limit T^Owe can easily see that 
the purely Gaussian results for AM and Cy calculated 
above are dominant. However, for the longitudinal sus- 
ceptibility the spin-wave interactions modify the purely 
Gaussian result. In this case, it is straightforward to show 
that for, H = H c2 , T -> 0, the dominant is x oc T 1 / 4 , 
instead of x oc T 1 / 2 calculated before. 
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FIG. 1: (Color online) Scaling plot in logarithmic scales of 
the magnetization data for the compound DTN obtained for 
fields up to 17T and temperatures T = 0.60,0.64,0.72 and 

0. 947^. The line shows the asymptotic behavior of the scaling 
function R(t) ~ t 3/2 in Eq. JT5), for H -> H c2 . The arrows 
indicate the region of validity of the scaling. 

Scaling analysis of the Magnetization — We start from 
the free energy density, which close to the zero tempera- 
ture quantum phase transition has the scaling form given 
by Eq. I|10p. The zero temperature critical exponents a, 
v and the dynamic exponent z are related to the dimen- 
sionality of the system d by the quantum hyperscaling re- 
lation, 2 — a = v(d+z)2&. In general for d e ff — d+z > 4, 

1. e., above the upper critical dimension d c = 4, the expo- 
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nents associated with the QCP at 5 = take Gaussian 
values, and in particular the correlation length exponent, 
v = 1/2. That this is the case in the present theory can 
be immediately verified writing the thermodynamic func- 
tions in a scaling form and identifying the relevant expo- 
nents. Furthermore Eq. (J6j) yields the dynamic exponent 
z = 2. Using the relation AM cx df/dH we get, 

AM oc\5\ 1 ~ a R{T/\S\» z ), (15) 

where AM = M sat (T, H sat ) - M (T, H) and M sat is mea- 
sured at the highest fields, H sat > 15T. Using the hyper- 
scaling relation for 3c? we obtain 1 — a = ^(3 + z) — 1 = 
(l/2)(3 + 2)- 1 = 3/2. 

Figure [1] shows the scaling plot of the magnetization 
for the compound DTN in fields up to 17T and for sev- 
eral temperatures. The magnetization data was obtained 
using a vibrating sample magnetometer adapted to be 
used in a 3 He cryostat. The external magnetic field, pro- 
duced by a superconducting coil, was aligned with the 
tetragonal axis of the sample, necessary condition to in- 
duce BEC. As shown in the figure the experimental data 
collapses in a good scaling plot when using the critical 
exponents appropriate for three dimensions. It can also 
be seen in Fig. Q]that for H — > H C 2, the scaling function 
R(t -> oo) ~ i 3/2 , such that, in this limit AM 3d cx T 3 / 2 
in agreement with the theory. The lower (upper) branch 
in Fig. Q] corresponds to data for higher (lower) fields 
than the critical magnetic field H C 2(0) =12.3T. This crit- 
ical field obtained from a criterion of best data collapse is 
in very good agreement with that found by Paduan-Filho 
et al.— using numerical differentiation of the magnetiza- 
tion data. We point out that a good scaling of the data 
is observed for magnetic fields sufficiently close to the 
critical field, i.e., for 11.92T< H c2 <12.6T. 

Summary — In spite of the large literature on the sub- 



ject of BEC of magnons, the phase diagram and the ther- 
modynamic properties around the upper critical mag- 
netic field H C 2, have not yet been completely examined. 
Since, in general, high magnetic fields are necessary to 
reach H C 2, experimental results and consequently theo- 
retical work are much more scarce in this region of the 
phase diagram. As we pointed out, the organic com- 
pound NiCh-ASC(NH 2 )2 (DTN) is ideal for this kind 
of studies since detailed magnetization curves can be ob- 
tained for very low temperatures close to H C 2- With this 
motivation we introduced a BOMFT approximation to 
study theoretically the upper field transition. We have 
obtained the dominant low temperature behavior of the 
magnetization (AM cx T 3 / 2 ), specific heat (Cy cx T 3 / 2 ) 
and susceptibility (x cx T 1 / 4 ) at the quantum critical 
trajectory and determined the shift exponent of the Neel 
line. We pointed out that, although the magnon-magnon 
interactions are irrelevant in the RG sense close to the 
QCP, they should be taken into account and determine 
the temperature dependence of the critical line and that 
of the susceptibility along the quantum critical trajec- 
tory. Our mean-field approach is justified since the effec- 
tive dimension for the transition at the QCP (H = H C 2, 
T = 0) is above the upper critical dimension. Finally us- 
ing the theoretical prediction we obtained for the scaling 
form of the field and temperature dependent magnetiza- 
tion close to H = H C 2, T = 0, we performed a scaling 
analysis of our magnetization data for DTN. The very 
good agreement between the theoretical and experimen- 
tal results provides unequivocal evidence that the tran- 
sition at H C 2 is a BEC of magnons. 
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